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ABSTRACT 


Results  of  an  experimental  Investigation  carried  out  in  the 
24-inch  water  tunnel  at  DTNSRDC  are  presented  for  the  cavltating 
propeller-induced  surface  pressures  and  a  reference-area  surface 
force  measured  on  a  flat  plate  boundary  representation  of  a 
nearby  ship  hull*  Newly  designed  apparatus  for  the  measurement 
of  surface  force  on  an  instrumented  disc  is  described  in  detail. 
It  involves  special  dynamometry  that  can  be  used  in  a  reverse 
mode  so  that  the  force  disc  can  be  driven  as  a  mechanical  shaker. 
Reciprocity  calibration  measurements  were  made  with  this  equip¬ 
ment  and  the  resulting  pressure-to-acceleration  transfer  function 
was  used  to  make  estimates  of  cavity  volume  velocity.  A  model  of 
a  seven-bladed  propeller  was  run  in  a  screen-generated  simulation 
of  the  steep  axial  wake  velocity  distribution  for  a  single  screw 
ship.  Results  are  displayed  for  the  longitudinal  and  transverse 
distribution  of  the  first  three  blade  rate  harmonic  components  of 
the  unsteady  pressure  amplitudes  and  the  unsteady  disc  force 
amplitudes  at  the  blade  tip  clearances  of  0.26D  and  0.4 ID,  over  a 
range  of  cavitation  numbers.  The  unsteady  pressure  pulse  ampli¬ 
tudes  obtained  in  this  work  compare  reasonably  well  with  measured 
values  obtained  previously  in  a  large  water  tunnel  experiment 
using  a  complete  scale  model  of  the  ship  hull.  The  experimen¬ 
tally  inferred  blade  rate  harmonic  component  of  cavity  volume 
velocity  from  the  reciprocity  measurement  agrees  approximately 
with  the  prediction  by  the  PUP-3  propeller  analysis  scheme. 


ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  Naval  Sea  Systems  Command  under  the  Ships,  Subs, 
and  Boats  Program,  Program  Element  62543N,  Task  Area  SF43421,  and  was  carried 
out  by  the  Ship  Performance  Department  under  Work  Units  1506-102  and  1506-103. 

INTRODUCTION 

Propeller-excitation  has  long  been  known  as  a  major  source  of  vibration  and 
noise  on  ships.  Problems  arise  in  the  form  of  excessive  hull  girder  vibrations 
in  the  stern  area  and  at  the  upper  levels  of  deckhouses;  unacceptable  localized 
vibrations  in  important  aft  end  spaces;  fatigue  damage  to  hull  plating,  internal 
stiffeners,  appendages,  or  other  structures  near  the  propeller;  high  levels  of 
Inboard  airborne  noise;  and  general  crew  nuisance.  Any  one  or  all  of  these 
symptoms  could  lead  to  the  imposition  of  speed  restrictions,  limitations  on  cer¬ 
tain  maneuvers,  or  avoidance  of  certain  propeller  RPM  regimes.  Understanding 
the  hydrodynamic  source  of  fluctuating  forces  and  unsteady  pressures  is  crucial 
for  the  eventual  prediction  of  the  magnitudes  and  spatial  extent  of  the 
excitation.  Thus  an  important  aspect  of  the  design  of  the  propeller-hull 


arrangement,  in  addition  to  satisfactory  steady  propulsive  performance,  is  the 
proper  accounting  for  the  propeller-induced  unsteady  forces,  moments,  and 
periodic  hull  pressure  fluctuations. 

There  are  two  main  categories  of  propeller  excitation.  Both  are  concerned 
with  the  interaction  of  the  propeller  blade  elements  operating  in  the  nonuniform 
velocity  wake  inflow  conditions.  (1)  Unsteady  shaft-transmitted  loads  (bearing 
1.  ads)  are  forces  and  moments  developed  on  the  propeller  and  transmitted  to  the 
hull  through  the  thrust  bearing  foundation  and  reduction  gearing  system  and  the 
aftermost  shaft  bearing.  The  most  important  class  of  problems  involving  unsteady 
bearing  loads  are  those  associated  with  the  longitudinal  shafting  vibration 
excited  by  the  unsteady  thrust  and  torque.  (2)  Unsteady  hull  surface  loads  or 
pressure-transmitted  loads  are  forces  and  moments  that  arise  from  the  spatially 
integrated  effect  of  t.ie  fluctuating  pressure  field  induced  by  the  propeller 
blades  passing  through  a  varying  inflow  velocity  pattern.  The  presence  of  inter¬ 
mittent  blade  cavitation  causes  further  magnified  pressure  fluctuations  due  to  a 
monopole  behavior  of  the  unsteady  cavity  volume.  Under  the  worst  conditions, 
periodic  pressure  pulses  having  nearly  constant  phase  angle  may  act  over  a  wide 
expanse  of  hull  surface  and  on  nearby  appendages.  In  this  situation,  integrated 
surface  forces  of  very  substantial  magnitude  can  be  delivered  to  the  hull. 

Presently,  methods  for  prediction  of  the  exciting  loads  are  different  for 
each  of  the  categories.  In  the  area  of  unsteady  bearing  loads,  prediction  capa¬ 
bility  has  been  available  for  some  time  and  It  has  been  well  integrated  into  the 
propeller  design  process  in  the  U.S.  Navy.  Notable  examples  of  analysis  and 
design  methods  and  verification  of  the  prediction  accuracy  are  provided  by 
Boswell  and  Cox**,  Valentine  and  Dashnaw^,  Boswell  ,  et  al^,  and  Boswell,  et  al^, 
and  in  the  numerous  fundamental  references  listed  in  these  works.  A  Z-bladed 
propeller  acts  as  effective  mechanical  filter  that  picks  out  only  certain  woke 
velocity  harmonic  components  and  their  radial  distributions  for  the  excitation  of  the 
various  unsteady  bearing  load  components  (Z  th  harmonic  for  fluctuating  thrust  and 
torque,  and  the  Z-l,  Z+l  harmonics  for  the  transverse  components  of  fluctuating  for¬ 
ces  and  moments).  The  application  of  bearing  load  analyses  to  design  has  been  rather 
successful  because  of  the  development  of  specific  design  modifications  such  as  blade 
skew,  warp,  etc  for  reducing  excessive  unsteady  blade  loads. 


*A  complete  listing  of  references  is  given  on  page  31. 
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Capability  for  the  analytical  prediction  of  unsteady  pressure  pulses  and 
distributed  surface  forces,  on  the  other  hand,  has  just  emerged.  There  is  a  vast 
and  growing  literature  on  this  subject  by  now,  largely  stemming  from  work  in  the 
early  1970’s  after  it  became  clear  that  unsteady  variations  of  cavity  volume  in 
conjunction  with  large  wake  ’’elocity  gradients  could  significantly  magnify  the 
blade  fanning  pressures  caused  by  a  propeller  operating  near  a  boundary.  A 
sampling  of  imoortant  references  concerned  with  prediction  of  unsteady  surface 
forces  and/or  periodic  hull  pressure  amplitudes  includes  Huse^,  Noordzij,® 

Vorus,  et  al?  Fitzsimmons,®  Kaplan  et  al^,  Hoshino,^  Kaplan,  et  al,^  Breslin, 
et  al*2f  and  Huse  and  Guoqiang.^ 

The  U.S.  Ncvy  has  been  fortunate  over  many  years  to  have  relatively  few 
experiences  with  the  severe  vibration  and  noise  problems  that  can  be  attributed 
to  surface  force  excitation.  Typically,  the  merchant  ship  designs  displaying 
problems  (tankers,  product  carriers,  Ro-Ro,  ferries,  LNG  ships)  have  tended  to 
exacerbate  the  difficulties  because  of  high  installed  power  (increasing  dramati¬ 
cally  through  the  lata  1960's  and  1970’s)  and  restrictive  propeller-hull 
arrangements,  giving  rise  to  steep  and  deep  wake  patterns.  In  contrast,  the 
open  stern  arrangements  of  most  Naval  combatants  and  many  auxiliary  ships  have 
wake  distributions  for  which  simple  vertical  tip  clearance  allowances,  such 
aa  the  values  recommended  by  Navy  customary  practice,  have  generally  been  satis¬ 
factory. 

Problems  of  excessive  interior  noise,  early  stage  blade  surface  erosion, 
and  heavy  localized  vibration  encountered  on  the  first  of  the  class  U.S.  Navy 
Auxiliary  Oiler  AO-177  have  been  described  by  Wilson,  et  al,^  and  illustrate  a 
case  where  the  neglect  of  surface  force  excitation  in  the  original  design  con¬ 
sideration  led  to  an  unsatisfactory  mismatch  of  propeller  and  hull  shape  (wake). 

In  this  case,  intermittent  blade  cavitation  resulted  in  excessive  unsteady 
pressure  pulses  over  a  wide  area  of  the  stern  near  the  propeller  and  unstable 
cavity  flow  behavior  near  the  blade  tips  and  trailing  edge  (cloud  cavitation). 

Among  the  many  lessons  learned  from  the  experimental  and  analytical  investiga¬ 
tions  carried  out  for  the  AO-177  project  was  the  observation  of  how  such  subtle 
changes  of  the  blade  cavitation  extent  and  appearance  could  produce  improvements 
in  the  excitation  magnitude.  Ar  important  result  of  the  improved  wake  achieved 
by  use  of  a  fin  on  the  AO  177  was  a  significant  reduction  of  magnitudes  of  unsteady 
pressure  pulse  amplitudes  distributed  nearby  and  away  from  the  propeller  plane. 
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The  hydrodynamics  of  the  blade  cavity  behavior  (thickness ,  total  cavity  volume 
variation,  and  interaction  with  tip  vortex  cavity)  appeared  to  be  crucial  to 
understanding  the  magnitudes  and  character  of  these  changes. 

Experiments  are  essential  for  providing  better  understanding  and  as  the 
source  of  empirical  estimates  of  realistic  excitation  magnitudes,  and  will  con¬ 
tinue  to  be  Important  because: 

a)  There  is  difficulty  and  uncertainty  in  predicting  analytically  detailed 
unsteady  blade  cavity  flows.  This  is  because  of  the  complex  three-dimensional 
tip  flow  around  propeller  blades  Involving  boundary  layer,  turbulent  and  tip 
vortex  flow  effects  coupled  with  the  unsteady  sheet  cavity  and  its  break-up  into 
a  bubbly  wake  flow. 

b)  There  is  a  need  to  predict  cavity  volume  velocity  at  all  Interesting 
harmonics,  not  only  at  the  blade  rate  value. 

c)  There  is  a  need  for  studying  other  Important  physical  effects  such  as 
the  Influence  of  wake  temporal  velocity  fluctuations  on  the  cavity  behavior  and 
coincident  propulsor-induced  pressures  and  forces  on  nearby  surfaces. 

The  purpose  of  the  present  project  was  to  design  and  build  an  experimental 
system  to  be  used  in  the  existing  water  tunnels  at  DTNSRDC  that  could  be  of  use 
in  obtaining  some  of  the  basic  information  desired  and  provide  a  possible  test¬ 
bed  for  design  assessment.  It  was  Intended  to  Include  enough  adaptability  to 
cover  standard  propeller-excitation  experiments,  as  well  as  new  investigations. 
The  set-up  was  also  directed  to  provide  some  measure  of  interim  capability  and 
development  of  experimental  skills  for  propeller-hull  Interaction  studies, 
before  the  advent  of  a  planned  new  large  facility^,  This  report  covers  the 
background  of  model  testing  on  propeller-excitation,  the  requirements  and 
details  of  the  experimental  set  up,  and  future  plans  for  the  system. 

BACKGROUND 

Model  experiments  are  very  Important  to  the  study  of  propeller-excitation 
characteristics  that  can  give  rise  to  ship  vibration  and  inboard  noise.  Many 
of  the  significant  milestones  in  our  understanding  of  the  phenomena  involved 
and  in  the  estimation  of  unsteady  loadings  for  full  scale  designs  of  propellers 
and  shafting  systems  have  been  achieved  through  model  testing  in  a  variety  of 


experimental  facilities.  It  is  anticipated  that  this  will  continue  to  be  true  as 
long  as  the  large  deficiencies  persist  in  the  reliability  (or  absence)  of  analy¬ 
tical  predictions  of  such  features  as  nominal  and  apparent  wake  velocity  pat¬ 
terns,  sheet  cavity  and  tip  vortex  cavity  geometry  and  volume  dynamics,  and 
interaction  of  hull-propulsor-wake-cavity  flows. 

Full  scale  experiments  are  also  essential  for  correlations  with  model 
results,  for  observing  flow  features,  and  for  crucial  data  on  important  measured 
ship  responses  such  as  localized  and  girder  vibration  and  inboard  noise. 

CHARACTERIZING  EXCITATION 

If  we  confine  our  attention  to  the  category  of  surface  pressure /surface 
force  and  moment  excitation,  we  can  characterize  the  propeller  exciting  magni¬ 
tudes  by  considering  several  main  features. 

1.  Fluctuating  pressure  amplitudes.  There  is  a  rotating  pressure  field 
associated  with  the  operation  of  a  propeller.  Figure  1  is  a  sketch  of  this 
pressure  disturbance  frozen  at  an  Instant  in  time  and  it  shows  a  definition  of 
the  spatial  phase  angle  measured  to  the  reference  line  of  the  blade  from  the 
surface  point  just  as  the  positive  peak  of  the  blade  rate  pressure  signal  passes 
(many  other  definitions  are  possible).  When  this  free  space  pressure  field 
interacts  with  a  rigid  nearby  boundary,  and  when  the  pressure  is  observed  from  a 
stationary  point  on  the  boundary,  the  amplitude  on  the  surface  will  fluctuate  in 
time  with  the  blade  passing  frequency  and  multiples  of  that  frequency.  The  sur¬ 
face  pressure  pulse  can  reach  up  to  twice  the  free  space  magnitude  because  of 
reflection  effect  from  the  boundary.  Point  pressure  measurements  on  a  nearby 
hull  surface  or  on  appendages  of  a  ship  or  a  model  can  be  obtained  experimen¬ 
tally  with  flush  mounted  pressure  transducers,  and  provide  some  picture  of  the 
pressure  loads  delivered  to  the  hull.  A  single  point  pressure  measurement  (over 
the  tip  of  the  propeller  at  the  closest  distance  to  the  hull)  is  sometimes  taken 
as  the  sole  indicator  of  the  magnitude  of  the  surface  force  excitation.  Such  a 
value  is  certainly  interesting  for  comparisons,  but  it  is  generally  insufficient 
to  characterize  a  representative  surface  force  magnitude  that  would  be  useful  to 
estimate,  for  instance,  the  hull  girder  vibration  response.  On  the  other  hand, 
the  distribution  of  point  pressures  over  a  region  of  hull  surface  may  well  be 
used  for  the  loading  excitation  for  structural  vibration  predictions. 
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In  any  case,  point  pressure  pulses  should  be  obtained  at  several  locations  for¬ 
ward  and  aft  of  the  propeller  plane  and  laterally  off  centerline  to  have  an  idea 
of  the  character  and  extent  of  the  loading. 

Some  example  of  measured  propeller-induced  pressure  pulse  distributions  are 
provided  by  Denny^  on  the  model  scale  and  by  Taniguchi  and  Ohtaka^  for  a  full 
scale  destroyer. 

2.  Reference-area  surface  force  amplitudes.  The  periodic  surface  force 
experienced  by  a  patch  of  area  of  the  nearby  hull  boundary  represents  a  net  or 
integrated  effect  of  the  distributed  pressure  pulses  exerted  on  the  surface. 

Such  a  quantity  can  be  obtained  experimentally  with  an  instrumented  surface 
segment  (See  for  example,  Kerwin,  et  al^)  or  by  integrating  measured  point 
pressures  over  the  desired  area,  accounting  for  the  spatial  and  temporal  (phase 
angle)  variations.  A  manageable  size  of  the  reference  area  might  be  on  the  order 
of  the  propeller  disc  area. 

3.  Surface  force  density.  The  longitudinal  distribution  of  the  induced 
surface  force-per-unit  length  is  called  force  density.  It  is  a  convenient  way 
to  represent  the  fluctuating  loading  imposed  on  the  hull  girder,  especially  in 
relation  to  excitation  of  overall  girder  vibration.  It  cannot  be  determined  by 
direct  measurement,  but  can  be  approximated  by  performing  girthwise  integration  of 
measured  point  pressures  resolved  as  vertical  or  other  force  components  at  each 
section  down  the  length  of  the  hull  shape.  It  is  noted  here  because  it  is  espe¬ 
cially  useful  for  displaying  the  rate  of  change  of  exciting  surface  force  as  a 
function  of  distance  in  the  vicinity  of  the  propeller,  and  away  from  it  as  well. 

4.  Net  surface  force  amplitude.  In  principle,  the  total  integrated  or  net 
surface  force  amplitudes  induced  by  the  propeller  can  be  determined  approxima¬ 
tely  from  area  integrations  of  measured  point  reassure  amplitudes.  There  has 
never  been  a  direct  model  experimental  measurement  made  of  the  total  induced 
vertical  surface  force  (separated  from  any  moment  effects),  and  none 
accomplished  on  a  full  scale  ship.  The  net  force  magnitude,  at  various  blade  rate 
harmonics,  would  represent  the  most  condensed  parameter  for  measuring  excitation 
level  important  to  girder  vibration,  although  of  less  importance  to  inboard 
noise.  The  net  vertical  force  is  not  likely  to  act  exactly  at  the  propeller 
plane  location  because  of  asymmetries  in  the  longitudinal  distribution  of  force 
density. 
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Example  discussions  of  vertical  surface  force  values  obtained  either  from 
distribution  of  measured  pressure  amplitudes  or  by  calculation  using  computed 
distributions  of  pressure  pulses  can  be  found  in  references  by  Hylarides,^ 

Skaar  and  Raestad^O,  and  Huse.^  Recommendations  on  the  levels  of  vertical 
force  ratio  (total  surface  force-to-mean  propeller  thrust)  that  could  cause  hull 
girder  vibration  problems  could  be  inferred  from  each  of  these  references,  but 
there  is  no  agreement  among  them.  Recommended  ratios  for  the  critical  value  of 
surface  force-to-thrust  ratio  for  the  onset  of  unacceptable  vibration  excitation 
seem  to  range  from  about  7  percent  to  30  percent.  Occasionally  values  well 
above  100  percent  have  been  predicted  with  no  report  of  excessive  ship  vibra¬ 
tion  on  the  ship  itself.  There  is  considerable  need  for  improvement  on  how  to 
interpret  and  use  estimated  values  of  net  surface  force  amplitude. 

5.  Cavity  volume  pulsation  strength.  For  a  propeller  operating  with  inter¬ 
mittent  blade  cavitation  in  a  nonuniform  flow,  the  resultant  pressure  amplitudes 
are  typically  increased  noticeably  over  noncavitating  levels.  The  increase  is 
relatively  even  greater  for  the  corresponding  amplitudes  of,  say,  the  vertical 
surface  force  density  or  the  total  surface  force  because  they  are  area- 
integrated  quantities  of  pressure  fluctuations  having  nearly  constant  phase 
angles  everywhere.  The  dominating  feature  is  the  monopole-like  pressure  source 
behavior  of  the  fluctuating  cavity  volume  of  each  blade  sheet  cavity.  Although 
the  primary  characteristic  of  a  typical  sheet  cavity  volume  variation  is  similar 
to  that  of  a  pulsating  spherical  bubble,^  there  are  co  itributions  from  cavi- 
tating  tip  vortices^*  and  possibly  significant  contributions  from  bulging  tip 
vortices. 22  The  term  cavity  volume  variation  refers  to  the  total  volume  within 
the  blade  cavities  plotted  as  a  function  of  the  circumferential  position  angle. 

A  useful  representation  of  the  main  physical  mechanism  involved  with 
exciting  pressures  from  cavity  pulsation  is  given  in  a  simple  expression  for  the 
pressure  fluctuation  produced  at  a  distance  r  from  the  center  of  a  pulsating 
spherical  bubble 
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where  the  cavity  volume  variation  is  given  by  ¥c(t).  This  shows  that  the  spa¬ 
tial  rate  of  decay  of  the  induced  pressure  amplitude  is  r”1  and  that  the  pulsation 

•  •  • 

strength  is  proportional  to  the  cavity  volume  acceleration  (or  for  harmo¬ 
nic  variations).  An  important  difficulty  in  making  use  of  this  type  of  rela¬ 
tionship  is  obtaining  accurate  determination  of  the  cavity  volume  variation, 
either  by  analytical  prediction  or  by  experiment. 

It  is  possible  to  determine  cavity  volume  pulsation  strengch  experimentally 
by  an  indirect  means  involving  the  concept  of  reciprocal  measurements.  From  the 
reciprocity  principle  for  linear  dynamic  systems,  a  simplified  expression  is 
available  (see  Appendix)  to  cover  the  case  of  the  unsteady  force  on  a  body  or  a 
body  surface  element  fixed  in  position  relative  to  a  fluctuating  point  source. 
Then,  in  a  water  tunnel  or  full  scale  test  environment,  the  "volume  velocity"  of  a 
fluctuating  blade  sheet  cavity  can  be  approximated  by 

4.  „  _IL _  (2) 

c  iw{p1/a2) 


where 


■V-c  -  time  rate  of  change  of  cavity  volume 

F2  *  amplitude  of  measured  unsteady  force  on  a  body  or  body  element, 
harmonic  component  corresponding  to  the  frequency  <*> 


m  -  circular  frequency 

(p^/a2)  ■  transfer  function,  measured  ratio  of  pressure  amplitude  at 

propeller  location  to  the  acceleration  amplitude  at  the  body 

i  »  (- 1)** 


All  the  quantities  in  this  formula  except  <*>  are  regarded  as  complex.  A 
discussion  of  this  expression  and  related  forms  is  given  in  the  Appendix. 

To  implement  the  reciprocity  technique,  there  are  two  phases  of  measure¬ 
ments.  During  the  calibration  phase,  the  body  element  is  forced  to  oscillate 
through  a  range  of  frequencies  of  interest.  The  transfer  function  ratio  (pi/a2) 
is  formed  for  each  frequency  from  the  simultaneous  measurement  of  the  pressure 
PI  measured  at  a  position  in  the  fluid  representing  the  approximate  center  of 
cavitation  and  the  vertical  acceleration  a2  of  the  forcing  body  element.  The 
final  measurement  phase  is  run  with  excitation  from  the  cavitating  propeller. 


8 


Then  once  the  cavitation-dependent  component  of  the  vertical  surface  force  blade 
rate  amplitude  (Fg)g  on  the  body  element  is  determined,  the  resulting  blade  rate 
component  of  volume  velocity  cap  be  evaluated  from  the  expression  given  above. 

REVIEW  OF  EXPERIMENTS 

It  is  useful  to  outline  the  scope  of  some  example  model  experimental 
investigations  that  have  been  used  to  explore  various  features  of  propeller- 
induced  surface  pressure/surface  force  excitation. 

Early  Experiments  on  Exciting  Loads  and  Induced  Pressures. 

The  earliest  work  that  showed  an  appreciation  of  the  presence  of  propeller- 
induced  surface  forces  in  addition  to  the  bearing  loads  was  that  of  Lewi8^3,24t 
An  experimental  scheme  was  devised  for  testing  rigid  ship  models  in  a  towing 
basin,  where  the  model  vibration  response  was  used  to  determine  the  exciting 
forces  under  self-propelled  conditions  by  comparing  the  response  to  that  pro¬ 
duced  by  a  calibration  vibration  exciter  located  in  the  stern  that  delivered 
known  oscillating  force  levels.  Lewis  was  able  to  decide  that  the  major  contri¬ 
bution  to  the  measured  propeller-induced  effective  vertical  force  on  the  hull, 
for  example,  was  due  to  the  surface  force.  Similar  measurements  were  made  on  a 
model  of  a  cargo  ship  form  and  reported  by  Lewis  and  Tachmindji .25  Continuing 
with  an  improved  rigid  model  approach,  but  with  smaller  models,  Lewis  presented 
a  summary  of  results  of  such  experiments  for  several  different  examples. 26 
Unsteady  blade  cavitation  played  no  role  in  these  model  basin  studies,  and  was 
not  suspected  as  being  important  until  many  years  later. 

An  experimental  program  to  measure  the  aft  end  vibratory  surface  loads  and 
inferred  bearing  loads  was  presented  by  Stuntz,  et  al^^.  In  this  case,  the 
stern  end  of  the  ship  model  was  mounted  on  flexures  and  instrumented  to  measure 
the  vertical,  horizontal,  and  torsional  unsteady  loads  induced  by  the  operating 
propeller  as  tested  in  a  towing  basin.  No  conditions  involving  intermittent 
blade  cavitation  could  be  considered. 

Early  experiments  conducted  in  water  to  measure  the  periodic  pressure 
disturbances  ahead  of,  behind,  and  leside  an  operating  propeller  have  been 
presented  in  reports  by  Tachmindji  and  Dickerson28,29  on  the  free  space 
pressure  pulses  near  a  propeller  running  in  uniform  flow  and  behind  a  strut;  by 
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Pohj30  on  the  fluctuating  pressures  Induced  on  a  flat  plate  surface  mounted 
parallel  to  the  propeller  axis  In  a  water  channel  and  on  the  pressure  pulses 
acting  on  the  hull  surface  of  a  model  ship  in  a  towing  basin;  by  Taniguchi3*  on 
the  pressures  Induced  on  a  flat  plate;  and  by  Weitendorf32  on  the  variations  of 
fluctuating  pressures  Induced  on  flat  and  concave  curved  plates  at  various  tip 
clearance  ratios,  on  pressure  pulses  on  flat  plates  with  various  elastic  moun¬ 
tings,  and  on  pressure  pulses  at  locations  near  the  propeller  on  the  hull  surface 
of  surface  ship  models.  All  these  early  tests  Involved  non-cavitating  flow 
conditions.  Experiments  by  Nelka33  were  also  concerned  with  the  Induced 
pressure  amplitudes  on  a  flat  plate  boundary  near  a  propeller  and  showed  the 
effects  of  tip  clearance  ratio  with  a  series  of  propellers  with  increasing  blade 
skew  angles  and  blade  warp  angles  and  the  effect  of  nonuniform  Inflow  velocity 
pattern. 

Examples  of  complete  ship  model  experimental  results  obtained  in  towing 
basin  tests  for  propeller-induced  vibratory  pressure  amplitudes  on  the  hull  sur¬ 
face  near  the  propeller  have  been  reported  by  Keil3^  for  an  oceanographic 
ship;  by  Huse35  for  a  tanker  form;  by  Weitendorf32  for  semi-tunnel  tunnel  stern 
combatant  form;  and  by  Jacobs,  et  al3*’  for  a  destroyer  hull  form.  All  these 
have  been  conducted  under  nonuniforra  inflows,  but  noncavitating  conditions,  so 
that  crucial  interaction  between  nonuniforra  inflow  and  intermittent  cavity 
volume  variation  is  missing  in  each  of  these  investigations. 

Experiments  with  Cavitating  Propellers  and  Modelled  Wake  Inflow 

Simple  Boundary  Arrangements.  Early  water  tunnel  experiments  were  carried 
out  by  Denny^,  who  made  extensive  measurements  of  the  propeller-induced 
periodic  pressure  amplitudes  on  a  flat  plate  boundary.  A  number  of  important 
effects  were  studied,  including  tip  clearance  ratio,  propeller  loading,  nonuniform 
inflow,  propeller  geometry  (expanded  area  ratio  and  blade  skew  angle),  and  cavi¬ 
tation  in  uniform  flow.  Extensive  comparisons  were  made  with  theoretical 
predictions  available  at  the  time,  and  these  were  used  to  check  the  separate 
contributions  of  loading  and  thickness.  It  is  interesting  that  Denny  was  the 
first  to  note  that  blade  cavitation  served  to  accentuate  the  magnitudes  of  the 
induced  pressure  pulses  in  the  case  of  uniform  inflow  and  extensive  sheet  cavi¬ 
tation.  The  increases  of  pressure  pulse  level  in  the  steady  cavitating  flow  were 
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comparable  to  the  Increases  caused  by  noncavltatlng  operation  in  certain 
of  the  wake  patterns  tested.  However,  the  crucial  combination  of  nonuniform 
flow  with  blade  sheet  cavitation  was  not  included  by  Denny,  so  the  important 
influence  of  intermittent  blade  cavitation  was  not  revealed  in  these 
experiments. 

The  pivotal  paper  by  Takahashi  and  Ueda®?  describes  measurements  of 
propeller-induced  pressure  amplitudes  on  a  flat  plate  in  both  noncavltatlng  and 
cavitating  conditions  and  in  both  uniform  and  nonuniform  wake  flow.  It  is  this 
work  that  is  often  credited  with  first  providing  the  experimental  motivation  of 
the  connection  between  intermittent  blade  cavitation  and  exaggerated  pressure 
pulses  that  arise  in  operation  in  nonuniform  flow. 

Other  water  tunnel  experiments  concerned  with  distributions  of  fluctuating 
pressure  amplitudes  induced  on  a  flat  plate  boundary  by  a  cavitating  model 
propeller  in  nonuniform  flow  have  been  carried  out  and  reported  by 
Weitendorf38»39  Takahashi, Chiba  and  Hoshino,^1  Sasajima,^*  Takekuraa,^  and 
Weitendorf^l  who  provides  a  summary  discussion  of  a  number  of  different  effects. 

Water  tunnel  experiments  involving  measurement  of  the  localized  propeller- 
induced  excitation  levels  on  somewhat  more  complicated  boundaries  have  been 
carried  out,  for  example,  by  Nelka^  on  the  pressure  pulses  produced  on  the 
inside  surface  of  a  duct  of  an  ducted  propeller  system,  and  by  Peck^^  0n  the 
pressure  pulses  at  points  on  the  centerline  of  a  tunnel  enclosing  65  percent  of 
the  diameter  of  two  different  propellers.  Kerwin,  et  al*®  presented  results  on 
the  levels  of  propeller-excited  unsteady  normal  surface  force  measured  on  a  flat 
force  disc  imbedded  flush  with  a  flat  boundary,  for  a  series  of  propellers. 

Partial  Stern  (Dummy  Model)  Experiments.  The  use  of  foreshortened  ship 
hull  models,  often  termed  dummy  models,  in  variable  pressure  water  tunnel 
experiments  with  operating  propellers  has  long  provided  a  practical  approach  to 
the  problem  of  attempting  to  model  the  proper  three  dimensional  wake  velocity 
patterns  for  model  propeller  cavitation  tests.  Since  the  early  1970's,  dummy 
models  have  also  been  used  in  many  laboratories  around  the  world  to  try  to  sur¬ 
vey  the  propeller-induced  periodic  pressures  at  points  near  the  propeller,  with 
approximate  representation  of  the  wake  velocities,  hull  boundary  shape,  and 
cavitation  number  environment. 
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Some  early  published  Investigations  using  foreshortened  hull  models  in  the 
determination  of  cavitating  propeller-induced  pressure  pulses  include  van  Manen,^^ 
Huse,^  van  Oossanen  and  van  der  Kooij,^®  and  van  der  Kooij  and  Jonk.^9 
Examples  of  more  recent  publications  that  have  either  made  use  of  results  from 
dummy  model  tests  or  mention  this  test  technique  in  some  fashion  are  Hylarides,5° 
Glover,  et  al,5*  Glover  and  Patience, 5 2  Lover  and  Wills, 53  Chiba,  et  al,  54 
Reed,  et  al,55  Lovik,5&  Munk,  et  al,57  and  Hadler,  et  al58. 

Useful  and  informative  work  on  the  test  techniques  and  interpretive  problems 
of  dummy  model  testing  in  water  tunnels  has  been  reported  on  by  English59  and  by 

Larsen. 60 

Complete  Model  Experiments.  Because  of  the  importance  of  the  wake  velocity 
pattern  to  the  correct  modelling  of  the  intermittent  blade  cavitation,  the  idea 
of  using  complete  ship  hull  models  for  cavitation  experiments  is  very  appealing. 
The  idea  is  that,  at  least  on  the  model  scale,  the  correct  wake  pattern  is 
represented  closely  by  using  the  geometrically  scaled  hull  form.  There  are  two 
prominent  facilities  actively  used  for  research  and  engineering  assessments  of 
ship  propulsion  systems  using  complete  ship  hull  models:  the  MARIN 
Depressurized  Towing  Tank  in  Wageningen  Netherlands,  which  provides  for  the 
simultaneous  satisfaction  of  the  Froude  number  and  the  cavitation  number;  and 
the  large  cavitation  tunnel  (Tunnel  No.  2)  at  SSPA  in  Sweden  which  provides  no 
free  surface  boundary  at  the  model  waterline,  but  employs  typical  free  stream 
velocities  that  lead  to  relatively  high  model  blade  Reynolds  numbers.  Good  sum¬ 
maries  of  some  experimental  results  and  discussions  of  testing  techniques,  data 
analyses,  and  comparison  studies  are  given  by  van  der  Kooij^l  on  the 
Depressurized  Towing  Tank,  by  Lindgren  and  Bjarne^  on  the  large  cavitation  tun¬ 
nel  at  SSPA,  and  by  Dyne  and  Hoekstra^3  0n  comparisons  of  results  between  the 
two. 


Reciprocity  Experiments.  Application  of  a  simplified  form  of  the  reciprocity 
principle  to  shipboard  measurement  of  the  propeller  source  strength  (volume  velo¬ 
city)  was  developed  in  References  64  and  65,  with  particular  reference  to  the 
problem  of  inboard  noise  transmitted  by  means  of  ship  structure  vibration 
excited  by  the  propeller.  The  basic  ideas  for  this  approach  are  discussed  in 
the  Appendix.  Gray66  has  also  demonstrated  the  use  of  the  reciprocity  concept 
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in  determining  the  source  strength  of  cavity  pulsation  and  the  estimate  of 
oscillating  hull  pressure  amplitudes  on  full  scale  ships.  Van  der  Kooij®!  has 
indicated  that  the  same  ideas  have  been  used  to  measure  the  model  scale  cavity 
source  strength  from  model  tests  in  the  NSMB  Depressurized  Towing  Basin. 

Vorus  in  a  discussion  of  a  ^aper  given  by  Kerwin,  et  al*®  proposed  ^n 
experimental /analytical  method  for  estimating  the  unsteady  force  on  a  body  (any 
shape)  located  near  a  fluctuating  point  source  by  exploiting  the  results  of  a 
reciprocity  experiment  conducted  with  very  simple  model  geometry  in  a  water  tun¬ 
nel.  This  suggestion  was  subsequently  followed  by  experiments  carried  out  with 
the  instrumented  force  disc  described  in  Reference  18,  using  a  model  of  a 
destroyer  propeller.  The  calibration  experiment  to  determine  the  acceleration- 
to-pressure  transfer  function  is  described  by  Whalen,  et  al.®^  Description  of 
the  wake  simulation  tests  and  results  of  the  measurement  of  the  disc  force 
amplitude  leading  to  the  determination  of  cavity  volume  velocity  are  presented 
by  Whalen,  et  al®®.  Application  of  these  results  to  determine  the  effect  of 
blade  tip  clearance  on  the  propeller-induced  vertical  hull  surface  force  for  the 
DO  963  has  been  summarized  by  Vorus. 

An  attempt  to  exploit  a  more  complete  statement  of  the  reciprocity  principle 
has  been  carried  out  by  des  Moutis^  in  the  M.I.T.  water  tunnel  for  the  experi¬ 
mental  determination  of  cavity  volume  velocity  as  well  as  the  six  unsteady 
propeller  bearing  force  and  moment  amplitudes. 


APPARATUS  AND  INSTRUMENTATION 


FACILITY 

The  particular  experimental  equipment  developed  for  this  work  was  designed 
for  use  in  the  24-inch  Variable  Pressure  Water  Tunnel  at  DTNSRDC.  This  facility 
provides  the  basic  propeller  dynamometry,  adjustable  ambient  tunnel  pressure  to 
simulate  cavitating  flow  conditions,  and  an  access  chamber  on  the  top  side  of 
the  open  jet  section  which  provides  space  for  the  force  measuring  apparatus. 


MOUNTING  PLATE  AND  INSTRUMENTATION 

A  one-inch  thick,  aluminum  flat  plate  insert  was  positioned  horizontally 
and  parallel  to  the  axis  of  the  open  jet  test  section  of  the  water  tunnel.  This 
functioned  as  the  representative  boundary  surface  for  the  initial  experiments. 
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With  the  propeller  installed  along  the  tunnel  centerline  on  the  downstream 
shaft,  the  vertical  position  of  the  plate  can  be  varied  from  7.5  in. (19.1  cm)  to 
9.0  in  (22.9  cm)  from  the  centerline  to  provide  vertical  tip  clearances  in  the 
range  of  az/D  =  0.25  to  0.4  for  a  nominal  10-inch  (25.4  cm)  diameter  model  pro¬ 
peller.  The  upstream  end  of  the  plate  was  faired  into  the  tunnel  convergence  with  a 
flexible  nylon  sheet  joining  the  plate  to  a  nose  piece  fastened  to  the  tunnel 
wall.  This  allowed  a  simple  adjustment  of  the  vertical  position  of  the  plate 
boundary. 

Force  Measuring  Apparatus 

The  surface  force  measuring  device  is  shown  in  the  sketch  of  Figure  2.  It 
consists  of  an  instrumented  force  disc  connected  to  a  robust  inner  plate  and 
backup  support  structure.  The  entire  support  system  is  connected  to  the  water 
tunnel  structure  through  the  flexures  of  two  block  gauges. 

The  force  disc  itself  is  a  circular  aluminum  plate  with  a  diameter  of  10  in. 
(25.4  cm)  and  thickness  of  1  in. (2. 54  cm)  mounted  flush  with  the  flat  plate 
boundary.  An  U-ring  was  fitted  in  the  gap  between  the  disc  and  the  hole  in  the 
mounting  plate.  The  disc  is  attached  to  the  inner  support  plate  by  a  tripod 
arrangement  of  three  piezoelectric  force  gauges.  The  particular  force  gauge 
(Wilcoxon  Model  L10)  was  chosen  because  it  satisfies  several  important  require¬ 
ments:  (a)  accurate  measurement  of  oscillatory  force  levels  with  amplitudes  expe¬ 

cted  in  the  range  of  5  x  10"2  to  5  lbf  (0.22  to  22  N)  in  a  frequency  range  possibly 
up  to  1500  Hz,  (b)  large  value  of  effective  spring  stiffness  constant,  and  (c) 
capability  of  being  driven  electrically  in  a  reverse  mode  in  order  to  be  used  as  a 
force  excitor  by  placing  an  oscillatory  voltage  signal  at  the  input. 

Various  masses  and  support  stiffnesses  of  the  elements  of  the  force 
measuring  device  were  chosen  to  obtain  a  flat  dynamic  response  characteristic  of 
the  force  disc  in  the  range  of  frequency  covering  the  first  three  harmonics  of 
the  blade  rate.  The  force  disc  mass  plus  its  added  mass  is  about  1  slug  (14.6 
kg)  mounted  on  the  force  gauges  having  a  combined  spring  stiffness  of  15  xlO^ 
lb/in.  (26.3  x  10^  N/cm).  The  massive  support  structure  includes  two  stacks  of 
weights  and  was  sized  to  provide  a  large  inertial  mass  of  about  550  lb  (17.1 
slug  or  249  kg)  suspended  on  the  flexures  having  an  effective  stiffness  of  1  x  10^ 
lb/in  (1.75  x  10^  N/cm).  The  system  was  designed  to  accommodate  the  exciting  f re- 
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quency  range  of  50  Hz  to  650Hz  which  brackets  the  range  of  the  first  three  blade 
rate  harmonics  for  any  model  propellers  of  interest ,  operating  at  the  necessary 
rates  of  rotation  (14  <n<  25  rps). 

The  force  disc  could  be  instrumented  with  as  many  as  three  accelerometers 
mounted  on  the  top  side  in  order  to  measure  the  vertical  vibration  levels  of  the 
disc  either  during  various  calibrations  or  during  an  excitation  experiment  with 
a  propeller  running  in  the  tunnel.  Waterproofing  of  the  accelerometers  and  all  the 
force  gauges  was  provided  by  coating  them  with  standard  potting  material.  Some 
details  of  Instrumentation  characteristics  are  given  in  Table  1. 

Surface  Pressures 

Mounting  holes  for  flush-mounted  pressure  gauges  were  arranged  in  the  pattern 
shown  in  Figure  3  in  the  flat  mounting  plate  and  the  force  disc.  There  are  52 
positions  on  the  flat  plate  and  13  on  the  force  disc.  These  provide  the  possi¬ 
bility  of  obtaining  a  comprehensive  picture  of  the  distribution  of  surface 
pressure  amplitudes  at  points  within  a  rectangle  of  dimension  1.1D  by  2D  in  the 
vicinity  of  the  propeller,  for  a  normal  10-inch  (25.4  cm)  diameter.  Specie! 
threaded  adaptor  fittings  permitted  the  use  of  either  of  the  two  different  types 
of  pressure  gauges  available  for  the  experiments.  The  expected  range  of  maximum 
fluctuating  pressure  amplitude  at  points  closest  to  the  tip  at  the  blade  rate 
frequency  or  any  of  its  higher  multiples  was  on  the  order  of  _+l  psi  (6.9  kPa). 

The  two  types  of  gauges  used  were  the  Kulite  XTM-190  and  the  CEC  4-312,  with 
nominal  operating  ranges  of  +10  psi  (69  kPa)  and  +25  psi  (172  kPa),  respec¬ 
tively.  Both  gauges  are  of  the  strain  gauged  diaphram  type.  Some  further  details 
on  the  pressure  gauges  are  given  in  Table  1 . 

CALIBRATIONS 

The  pressure  gauges  were  calibrated  in  a  static  pressure  test  stand  to  deter¬ 
mine  the  sensitivity  (slope  *  pressure/volt)  for  each  gauge.  These  slopes  were 
also  checked  at  intervals  throughout  the  experimental  program,  using  tunnel 
pressure. 

Extensive  calibrations  were  performed  on  the  surface  force  measuring  device 
in  order  to  establish  its  characteristics  as  an  unsteady  dynamometer  and  for  its 
application  in  a  reciprocity  experiment. 
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It  was  necessary  to  establish  the  correct  force  slope  factors  for  the  three 
force  gauges  and  to  determine  the  frequency  response  characteristic  for  the  force 
disc*  An  unsteady  force  was  applied  to  the  disc  by  means  of  an  electrodynamic 
shaker.  The  force  applied  to  the  disc  was  measured  independently  with  a  pre¬ 
calibrated  piezoelectric  force  gauge.  For  the  case  of  the  tunnel  full  of  water  and 
with  the  0-ring  in  place  around  the  edge  of  the  disc,  the  measured  dynamic 
response  is  shown  in  Figure  4,  displaying  the  ratio  force  output-to-f orce  input 
plotted  versus  frequency  of  excitation.  The  frequency  response  curve  is  essen¬ 
tially  flat  to  around  650  Hz  where  it  begins  to  rise  toward  a  natural  resonance 
frequency  at  930  Hz.  This  indicates  that  the  desired  level  response  characteristic 
holds  satisfactorily  out  to  frequencies  well  beyond  the  third  blade  rate  harmonic 
for  model  propellers  with  as  many  as  seven  blades.  The  flat  force  response  also 
permitted  the  final  determination  of  the  calibration  slopes  for  the  three 
piezoelectric  force  gauges  with  the  fully  assembled  dynamometer. 

The  second  major  function  of  the  force  disc  apparatus  is  to  provide  the 
possibility  of  determining  the  characteristic  cavity  volume  variation  of  a  cavi- 
tatlng  propeller  in  nonuniform  inflow  by  means  of  the  indirect  method  exploiting 
the  reciprocity  principle.  To  make  this  type  of  measurement  it  is  necessary 
first  to  perform  a  reciprocity  calibration  which  determines  the  transfer  function 
relating  the  fluid  pressure  response  at  a  field  point  to  the  magnitjd-*  of  ver¬ 
tical  vibration  of  the  surface  area  element  (force  disc).  The  reciprocity  calibra¬ 
tion  was  carried  out  by  oscillating  the  force  disc  vertically  along  its  axis 
(normal  to  the  surface)  by  means  of  a  sinusoidally  varying  voltage  placed  across 
the  three  piezoelectric  force  gauges  which  support  the  disc.  This  effectively 
drove  the  force  gauges  in  a  reverse  mode  so  that  the  disc  was  itself  the 
electrodynamic  shaker  for  the  calibration  phase.  Excitation  frequency  covered  a 
wide  range  from  about  70  Hz  to  1000  Hz  in  order  to  include  the  frequency 
variation  up  to  the  third  harmonic  of  blade  rate.  The  vertical  acceleration  of 
the  disc  was  measured  with  accelerometers  mounted  on  the  back  side  of  the  disc. 
Pressure  response  in  the  fluid  was  measured  in  the  absence  of  the  propeller 
using  a  hydrophone  located  at  a  point  in  the  12  o'clock  position  at  a  distance 
0,9  R  from  the  axis.  This  location  corresponds  roughly  to  the  expected  center 
of  the  largest  blade  cavity  volume  determined  from  visual  observation. 

The  resulting  calibrated  transfer  function  is  shown  in  Figure  5  as  the  curve 
of  the  logarithm  of  the  ratio  of  pressure  araplitude-to-vertical  acceleration 
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(p/a)  plotted  versus  frequency.  There  18  a  variation  with  frequency  apparent 
in  this  empirical  function,  with  a  local  peak  occurring  at  around  500  Hz,  and  a 
relative  flattening  beyond  800  Hz. 

For  frequencies  below  300  Hz  the  measured  values  of  the  transfer  function 
for  (p/a)  were  found  to  display  a  temporal  variation.  Since  the  reciprocity 
calibration  was  obtained  from  a  long  analog  tape  record  of  the  simultaneous 
measurement  of  fluid  pressure  and  disc  acceleration  taken  during  the  same 
calibration  session,  several  values  of  (p/a)  at  a  given  frequency  could  be 
calculated  from  several  different  segments  of  the  tape.  For  the  data  displayed, 
each  segment  consisted  of  a  frequency  sweep  through  the  range  of  interest.  This 
provides  a  check  on  repeatability  and  shows  that  the  values  of  (p/a)  were  reliably 
repeatable  for  the  frequency  regime  higher  than  300  Hz.  Below  300  Hz,  however, 
there  was  a  range  of  values  computed  for  (p/a)  at  any  one  frequency,  indicating 
that  the  transfer  function  varied  with  time.  In  Figure  5,  the  shaded  region 
shows  the  range  of  values  of  the  transfer  function  computed  from  several  dif¬ 
ferent  segments  of  the  measured  record. 

At  this  time,  there  is  no  conclusive  explanation  for  the  variability  of  the 
measured  transfer  function  in  the  low  frequency  range.  One  possibility  is  that 
under  unsteady  excitation,  there  are  pressure  reflections  from  structural  ele¬ 
ments  Inside  the  tunnel  and  test  section,  such  as  the  nozzles  or  the  installed 
flat  plate.  There  could  be  lnterferring  unsteady  pressures  set  up  by  induced 
vibrations  of  the  tunnel  walls  and  model  structure.  A  likely  contributor  to  the 
problem  could  be  poor  signal-to-noise  ratio  for  the  pressure  response  to  the 
disc  oscillation.  This  is  expected  to  be  worse  at  the  low  frequencies  because 
of  the  very  small  acceleration  levels  being  produced  by  the  mechanical  system  in 
the  frequency  range  involved.  It  is  not  likely  that  the  variability  of  the 
(p/a)  ratio  is  attributable  to  changing  local  properties  such  as  temperature  or 
ambient  pressure,  because  the  time  scales  of  these  changes  are  much  longer  than 
characteristic  time  involved  with  the  changes  of  the  transfer  function  magni¬ 
tude.  This  problem  requires  further  study. 

For  the  model  propeller  experiments  reported  on  here,  the  blade  rate  fre¬ 
quency  at  the  convenient  operating  speed  was  98  Hz.  At  this  frequency,  the 
variability  of  the  (p/a)  involves  a  range  of  about  +34  percent. 
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EXCITATION  EXPERIMENTS  WITH  MODEL  PROPELLER 

PROPELLER 

For  comparison  purposes,  it  was  decided  to  use  the  9.812-in.  (24.92  cm) 
diameter  DTNSRDC  model  propeller  number  4677  which  represents  the  21-ft,  seven- 
bladed,  skewed  propeller  for  the  single  screw  Naval  Auxiliary  Oiler  AO-177  with 
the  scale  ratio  25.682.  The  wake  adapted  design  of  this  propeller  is  described  by 
Valentine  and  Chase. Table  2  is  a  summary  of  the  particulars  of  the  propeller 
geometry.  Drawings  of  the  blade  shape  are  given  in  Figure  6.  Figure  7  shows 
the  propeller  open  water  characteristics  presented  by  Hendrican  and  Remmers^. 
Appropriate  operating  conditions  for  the  propeller  were  selected  on  the  basis  of 
the  model  powering  experiments  reported  in  Reference  72. 

EXPERIMENTAL  PROCEDURE 
Wake  Simulation 

All  the  experiments  shown  here  were  carried  out  in  a  non-uniform  flow  in  the 
tunnel  produced  by  a  wake  screen  which  provides  for  control  of  the  axial  velo¬ 
city  component  only.  The  wake  screen  was  designed  to  simulate  the  nominal  wake 
distribution  Vy/V  produced  by  the  AO-177  as  determined  from  a  towing  basin  wake 
survey  with  a  scaled  model  of  the  ship  hull.  Wake  surveys  were  conducted  in  the 
tunnel  with  the  screen-generated  flow  using  a  rake  of  five,  5-hole  spherically- 
headed  pitot  tubes.  The  flat  plate  boundary  was  set  at  two  locations 
corresponding  to  propeller  tip  clearances  of  az/D  ■  0.26  and  0.41.  Figures  8, 

9,  and  10  are  example  comparisons  of  the  measured  axial  velocity  ratio  Vy/V  ver¬ 
sus  9  for  r/R  ■  0.557,  0.774,  and  1.178  respectively,  showing  data  obtained  from 
a  towing  basin  wake  survey^*^  and  the  tunnel  wake  surveys.  Definition  of  the 
position  angle  9  is  given  in  the  figures.  These  comparisons  show  that  the  Vy/V 
flow  patterns  behind  the  wake  screen  for  the  two  different  plate  locations  are 
generally  similar  to  the  target  nominal  velocity  distributions,  but  are  not 
Identical.  One  noticeable  feature  is  that  there  are  differences  in  the  level  of 
the  Vy/V  values  when  the  position  angle  9  is  outside  the  main  velocity  defect 
region  (which  is  contained  inside  -90  deg  <  9  <  90  deg).  The  relatively  higher 
values  for  the  case  az/D  =  0.26  are  attributable  to  relatively  increased  average 
flow  velocity  caused  by  the  smaller  flow  cross  section  area.  This  speed-up 
influence  is  less  pronounced  for  the  case  of  az/D  =  0.41,  because  the  flow  area 
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blockage  Is  less.  The  narrow  velocity  defect  from  the  lower  skeg  region,  cen¬ 
tered  about  the  position  angle  9  ■  180  deg,  and  seen  clearly  In  the  towing  basin 
velocity  distribution,  was  not  reproduced  in  the  wake  screen  velocity  pattern. 

The  skeg  wake  velocity  feature  In  this  case  is  much  less  important  than  the  main 
hull  wake  characteristic  because  there  Is  no  known  blade  cavitation  that  occurs 
In  that  region.  For  this  reason,  dense  screen  material  was  not  applied  along  the 
centerline  of  the  lower  disc  area  of  the  wake  screen  used  for  the  present 
experiments. 

There  are  important  similarities  between  the  wake  screen  velocity  patterns 
and  the  towing  basin  velocity  distributions  (the  target  wake).  The  values  of 
the  minimum  V^/V  ratio  in  the  vicinity  of  9  ■  0  are  nearly  the  same  in  the  two 
wake  screen  distributions  as  in  the  target  wake,  although  the  precise  angular 
positions  may  be  offset  somewhat  for  the  wake  screen  cases.  More  importantly, 
the  relative  slopes  of  the  main  velocity  defect  with  respect  to  the  circumferen¬ 
tial  angle,  denoted  by  d(Vx/V)/d9,  are  approximately  the  same  for  the  target 
distribution  and  the  wake  screen  distributions  along  the  rays,  say  9  *  40  deg 
and  9  *  320  deg.  These  similarities  are  concerned  with  features  of  the  wake  that 
are  most  important  to  the  production  of  unsteady  pressures  due  to  cavity  volume 
variation,  as  noted,  for  instance  in  Reference  75.  It  is  concluded  that  pro¬ 
peller  blade  cavities  should  reach  approximately  the  same  3ize  at  about  the  same 
circumferential  position,  and  should  grow  and  diminish  at  about  the  same  rate  in 
either  of  the  two  wake  screen  velocity  patterns  as  in  the  towing  basin  velocity 
distribution.  Since  the  dominant  contribution  to  the  propeller  unsteady 
pressures  and  surface  force  comes  from  the  unsteady  cavity  volume  variation, 
experiments  in  the  wake  screen  velocity  patterns  described  above  could  be 
expected  to  provide  excitation  characteristics  generally  meaningful  to  the  case 
of  the  AO-177  propeller  and  wake  combination. 

Unsteady  Pressures  and  Disc  Surface  Force 

Measurements  have  been  carried  out  in  the  24-inch  water  tunnel  using  the  newly 
developed  apparatus  and  propeller  Model  4677.  The  propeller  was  mounted  on  the 
downstream  shait,  with  the  boundary  plate  arranged  parallel  to  the  axis. 

Pressure  gauge  measurements  were  taken  at  11  of  the  locations  on  the  disc  and  at 
many  locations  on  the  surrounding  boundary  plane,  concentrating  on  the  points 
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along  the  longitudinal  and  transverse  centerlines  (xpg  and  y  axes,  respectively) 
as  depicted  in  Figure  3.  The  propeller  reference  center  identified  in  Figure  6 
was  placed  at  the  longitudinal  position  of  the  transverse  axis,  xp^  =*  0.  Thus, 
the  location  of  a  plane  through  the  blade  tips  was  at  a  distance  0,1 15D 
downstream  of  the  xpg  »  0  position. 

Unsteady  disc  force  measurements  were  made  with  the  three  piezoelectric 
force  gauges  in  place,  and  in  some  cases,  simultaneously  with  pressure  gauge 
measurements. 

All  the  experiments  conducted  for  this  work  were  performed  at  the  loading 
condition  corresponding  to  the  advance  coefficient  J  *  0.79  and  thrust  coef¬ 
ficient  K-y  -  0.295,  producing  a  thrust  coefficient  of  Cfjj  -  1.2.  Test  section 
velocity  was  set  for  the  target  J  by  using  the  thrust  identity  method  from  the  open- 
water  characteristics.  This  was  accomplished  by  setting  the  shaft  rotation  speed 
for  the  desired  advance  coefficient  and  adjusting  the  tunnel  water  speed  until 
the  thrust  coefficient  matched  the  open-water  thrust  coefficient.  The  velocity 
calibrated  in  this  way  was  held  constant  and  the  tunnel  pressure  varied  to  change 
the  cavitation  number.  Shaft  rotation  speed  was  maintained  at  14  rps  for  all  the 
conditions  tested.  The  Reynolds  number  based  on  chordlength  at  the  0.7R  and  the 
approximate  total  velocity  was  0.51  x  10^. 

The  range  of  cavitation  numbers  included  on  *  17,86  to  8.36  (°va  "  28.6  to 
13.4)  which  covered  non-cavitating  reference  conditions  as  well  as  conditions 
with  mixed  tip  vortex  and  sheet  cavitation  characteristic  for  this  propeller.  The 
minimum  cavitation  number  used  here  was  limited  by  a  tunnel  flow  speed  that  would 
not  damage  the  wake  screen.  With  the  loading  conditions  held  the  same 
throughout  the  range  of  on,  the  effects  of  different  blade  cavitation  patterns  on 
the  unsteady  pressures  and  disc  force  could  be  studied  independent  of  the 
propeller  loading  condition. 

Total  air  content  of  the  tunnel  water,  measured  with  a  Van  Slyke  apparatus, 
was  held  in  the  range  of  24  to  35  percent  of  saturation  at  atmospheric  pressure, 
mainly  to  provide  clear  visibility.  Water  temperature  varied  from  69°  to  74°  F 
(20.60  t0  23.3°  C). 

Measurements  of  both  the  pressure  fluctuations  and  disc  unsteady  force  were 
obtained  with  the  flat  plate  boundary  positioned  at  two  different  locations 
corresponding  to  propeller  tip  clearance  ratios  az/D  =  0,26  and  0.41,  The  ver¬ 
tical  tip  clearance  ratio  on  the  AO-17?  ship  is  az/D  =  0.292. 
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’’’he  test  procedure  consisted  of  recording  the  pressure  and  force  transducer 
signals  for  several  minutes  at  a  particular  cavitation  number,  with  the  simulta¬ 
neous  recording  of  the  data  onto  magnetic  tape.  Throughout  most  of  the  test 
conditions,  the  pressure  and  force  signals  were  cheeked  online  with  a  time 
series  analyzer  to  review  the  character  of  the  pressue  and  force  spectra  and  the 
relative  contributions  at  the  blade  rate  frequency  and  its  higher  harmonics. 

Taped  pressure  and  force  gauge  records  consisted  of  a  recorded  digitized  signal 
at  256  points  around  each  shaft  rotation  as  triggered  from  a  shaft  encoder 
attached  to  the  prooeller  shaft.  A  timing  channel,  recorded  simultaneously  with 
the  pressure  and  force  signals,  provided  a  single  pulse  per  shaft  rotation  that 
was  arranged  to  indicate  the  passage  of  a  blade  reference  line  past  the  12 
o'clock  position  (point  of  closest  approach  to  the  flat  boundary).  It  was  used 
for  the  measured  phase  analysis. 

After  the  recording  sessions,  the  data  were  analyzed  using  a  Model  70 
Interdata  computer  to  provide  the  Fourier  Series  1  armonic  components  over  a  con¬ 
venient  specified  time  interval,  usually  10  seconc’s.  Thus  a  typical  analysis 
record  for  these  experiments  covered  140  propeller  revolutions. 

It  was  clear  from  the  harmonic  analysis  and  f  -'Mn  the  time  series  analyzer 
spectra  that  the  signals  of  the  surface  pressures  and  disc  force  were  dominated 
by  the  blade  rate  frequency.  Amplitudes  of  higher  multiples  of  blade  rate  were 
always  smaller  than  the  blade  rate  itself,  but  not  necessarily  with  magnitude 
decreasing  strictly  with  increasing  harmonic  order. 

Data  Analysis 

The  unsteady  behavior  of  the  pressure  measured  at  any  location  is  a  periodic 
fluctuation  p(t)  and  can  be  represented  by  the  Fourier  series. 


CO 

p(t)  =  y  (  A  cos  mcot  +  B  sin  mwt  )  (3) 

6 ,  m  m 

m=l 


rewritten  in  a  form  to  display  total  amplitude  and  phase  angle  as 
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p(t)  -  I  Pm  sin  («t>  -  Ym) 
m=l 


where  ♦  ■  is  the  position  angle  taken  positive  clockwise  from  12  o'clock, 
looking  upstream 

u  .  2*n  *  shaft  circular  frequency 

Pm  "  +  Bm^ )  */2  ■  total  amplitude  of  the  mth  harmonic  component  of  the 

pressure  fluctuation 

Ym  *  tan”*  (-Ajj/Bjj)  -  harmonic  phase  angle  of  the  mth  pressure  harmonic 

The  positive  peak  (maximun  positive  value)  of  the  mth  harmonic  component  of 
pressure  fluctuation  occurs  when 

"  Y  *  ~  (5) 

m  2 


Thus  the  spatial  phase  angle  at  which,  say,  the  Zth  or  blade  rate  harmonic  of 
pressure  reaches  a  maximum  positive  value  is 
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This  phase  angle  Is  the  angular  location  of  the  blade  reference  line  when  the 
blade  rate  pressure  positive  peak  occurs,  measured  clockwise  looking  upstream 
(see  Figure  1). 

The  amplitude  of  the  mth  harmonic  pressure  component  (Ap)m  *  Pm  is  made  non- 
dimensional  in  the  following  definition 
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In  the  present  work,  we  confine  our  Interest  to  the  results  for  the  first  three 
blade  rate  harmonics  ,  with  m  =  yZ 

m  *  Z,  2Z,  3Z 
so  that  y  =  1,  2,  3 


Similarly,  the  unsteady  reference  area  surface  force  on  the  disc  can  be 
expressed  in  terms  of  its  Fourier  series 


Fq(t)  *  I  (F  )  sin  (tn<j>  -  B  ) 

S  m=l  S  m  m  (8) 


where  4>  =  tot  is  the  position  angle 

(F  )m  *  total  amplitude  of  the  mth  harmonic  component  of 

fluctuating  disc  force  (obtained  from  the  sum  of  the 
three  force  gauge  signals) 

B  *  harmonic  phase  angle  of  the  mth  disc  force  component 
m 

The  nondimensional  form  of  the  disc  force  amplitude  for  the  mth  harmonic  is 
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Using  this,  the  ratio  of  the  unsteady  force  amplitude  to  propeller  mean  thrust, 
say  at  blade  rate  frequency,  is  the  ratio  of  the  coefficients 

_  (IVz 

T  Kt  (10) 


RESULTS 

Propeller  Cavitation 

Visual  observation  of  the  propeller  operated  under  the  test  conditions 
described  earlier  showed  that  tip  vortex  cavitation  springing  from  the  blade  tips 
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appeared  at  around  <?n  -  13.  As  the  cavitation  number  was  decreased  from  this 
value,  the  tip  vortex  increased  in  size.  Sheet  cavitation  began  forming  on  the 
suction  side  near  the  blade  tips  at  0n  near  12.  At  on  ■  10.5  there  was  substan¬ 
tial  sheet  cavitation  on  each  blade.  For  the  cavitation  number  <*n  *  8.36,  the 
lowest  value  achieved  in  this  test  series,  the  sheet  cavity  on  each  blade 
covered  a  portion  of  the  outer  25  percent  of  the  propeller  radius.  The  circum¬ 
ferential  extent  of  the  blade  sheet  cavitation  was  from  about  90°  before  12 
o'clock  around  to  60°  beyond  the  top-dead-center  position. 

Unsteady  Pressures 

Figures  11  and  12  show  the  variation  of  pressure  coefficient  values  Kp  versus 
cavitation  number  for  the  amplitudes  of  fluctuating  pressure  acting  at  the  center 
of  the  disc  (directly  over  the  propeller  center)  for  the  first  three  harmonics  of 
blade  rate.  Figure  11  shows  data  for  the  tip  clearance  ratio  az/D  *  0.26,  with 
Figure  12  applicable  to  the  larger  tip  clearance  az/D  -  0.41.  The  pressure  ampli¬ 
tude  curves  are  flat  in  the  nc-n-cavitating  regime  of  large  on.  p0r  values  of  c»n 
<  12,  the  pressure  amplitude  characteristics  start  to  rise  with  decreasing  cavi¬ 
tation  number.  The  noticeable  Increase  in  the  levels  of  all  the  (Kp)m  values 
coincides  with  the  appearance  and  growth  in  size  of  blade  sheet  cavities,  a  fact 
that  confirms  our  general  understanding  of  the  importance  of  unsteady  cavity 
volume  variation  to  the  pressure  excitation. 

For  both  tip  clearances,  the  blade  rate  harmonic  component  is  typically 
larger  than  the  second  and  third  blade  rate  contributions.  The  differences  are 
clearl'  dependent  on  the  tip  clearance  ratio.  At  the  smallest  tip  clearance  of 
az/D  1  0.26  the  blade  rate  amplitudes  exceed  the  higher  harmonics  by  a  factor  of 
more  than  2.5.  For  the  larger  az/D  *  0.41,  the  higher  harmonic  amplitude  values 
are  ’.oraparable  to  the  fundamental  blade  ra.e  amplitude. 

In  a  comparison  of  the  results  of  magnitudes  of  the  blade  rate  component  for 
the  wo  tip  clearance  ratios,  there  appears  to  be  a  large  influence  of  tip 
cleat ance.  The  smaller  tip  clearance  produces  the  larger  (Kp)z  levels  in  both 
the  ncn-cavitating  and  cavitating  regimes.  This  is  certainly  consistent  with 
previous  experimental  findings.  Some  of  the  difference  i:«  pressure  amplitudes 
between  the  two  tip  clearance  cases  may  be  attributed  to  differences  in  the 
screen-generated  nominal  wakes  for  the  two  tip  clearance  ratios,  as  seen  in 


24 


Figures  8  through  10.  The  extremes  of  the  hydrodynamic  advance  angle  away  from 
the  mean  value  are  different  for  the  two  wakes ,  and  this  could  influence 
the  variation  of  blade  cavitation.  However,  visual  observation  of  the  blade 
cavitation  for  the  two  tip  clearance  cases  revealed  no  significant  differences. 

It  is  plausible  that  for  the  gauge  location  over  the  propeller  center,  the 
measured  differences  of  pressure  amplitude  are  due  mainly  to  the  tip  clearance 
variation. 

The  longitudinal  distributions  of  the  first  three  blade  rate  harmonics  of 
fluctuating  pressure  amplitude  coefficient  along  the  centerline  of  the  plate 
(the  xpc“axis  shown  in  Figure  3)  are  displayed,  respectively,  in  Figures  13,  14, 
and  15  for  tip  clearance  ratio  az/D  ■  0.26.  The  comparable  distributions  of 
(K„)2  versus  xp^/D  for  the  tip  clearance  case  az/D  -  0.41  are  given  in  Figure  16 
through  18.  Here,  positive  values  of  xpg/D  are  upstream  of  the  propeller 
center.  The  contours  in  these  groups  of  figures  are  for  five  different  values  of 
cavitation  number.  These  plots  show  the  expected  fall-off  of  pressure  amplitude 
both  upstream  and  downstream  of  the  location  of  the  propeller  center. 

The  transverse  distributions  of  the  first  three  blade  rate  harmonics  of 
pressure  coefficient  are  plotted  versus  the  nondimens ional  distance  y/D  (along  the 
lateral  axis  in  line  with  propeller  center)  in  Figures  19  through  21  for  az/D  - 
0.26;  and  in  Figures  22  through  24  for  az/D  ■  0.41.  Here  the  positive  values  of 
y/D  are  taken  to  starboard,  looking  upstream.  Generally,  the  largest  blade  rate 
pressure  amplitudes  occur  near  the  propeller  centerline  and  fall  off  to  either 
side.  It  has  been  observed  from  published  experimental  work  with  a  variety  of 
propeller  types,  that  the  transverse  distributions  of  blade  rate  pressure  amplitu¬ 
des  are  not  necessarily  expected  to  be  symmetric  about  the  propeller  centerline 
or  about  a  line  passing  through  the  location  of  the  peak  pressure.  Propeller 
blade  skew  will  tend  to  accentuate  the  non-symmetry,  particularly  for  the  lower 
range  of  an  where  transient  cavitation  begins  to  dominate  excitation  pressures. 

The  propeller  used  in  the  present  experiments  has  45  deg  skew  at  the  tip  radius. 
With  blades  having  positive  skew  ("skew  back"),  even  with  a  symmetrical  wake, 
the  region  of  greatest  cavitation  extent  is  displaced  around  to  the  right  hand 
side  of  the  disc  area  for  a  right-handed  propeller  (see  for  instance,  Reference 
14).  For  very  low  cavitation  numbers,  the  largest  unsteady  pressures  are  asso¬ 
ciated  with  the  cavity-collapse  portion  of  the  cavity  volume  variation,  which 
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occurs  over  in  the  right-hand  quadrant  of  the  propeller  disc  rather  than  at  the 
12  o'clock  position  where  the  tip  clearance  is  smallest.  The  displacement  and 
distortion  of  cavity  volume  due  to  effects  of  skew  can  influence  the  higher  har¬ 
monics  of  the  surface  pressure  in  non-trivial  ways  as  well.  In  the  present 
experiments,  carried  out  at  intermediate  values  of  cavitation  number,  the  loca¬ 
tion  of  the  peak  blade  rate  pressure  amplitude  is  apparently  influenced  most  by 
the  effect  of  tip  proximity,  and  occurs  at  or  near  the  position  of  closest 
approach  of  the  blade  tip. 

For  higher  harmonics  of  blade  rate,  the  curves  of  pressure  amplitudes  at 
constant  on  -  values  are  sometimes  skewed  to  one  side,  or  may  show  some  slightly 
mixed  trend  with  respect  to  cavitation  number  on  one  side  of  the  centerline  that 
differs  from  the  opposite  side.  (See,  for  instance,  Figures  17,  20,  21,  23,  and 
24).  These  are  the  manifestations  of  a  combination  of  higher  harmonic  distortion 
by  a  skewed  cavitation  extent  (noted  above)  and  scatter  and  inaccuracy  asso¬ 
ciated  with  the  relatively  low  levels  of  pressure  magnitude  of  the  higher  har¬ 
monic  components  compared  with  the  blade  rate  contributions. 

Unsteady  Disc  Force 

The  measured  fluctuating  disc  force  amplitude  coefficients  for  the  first 
three  blade  rate  harmonics  are  shown  plotted  versus  cavitation  number  for  az/D  ■ 
0.26  in  Figure  25,  and  for  az/D  *  0.41  in  Figure  26.  The  symbols  plotted  in 
these  graphs  represent  the  average  values  for  two  or  three  data  spots  at  each  of 
the  tested  cavitation  numbers.  Similar  to  the  trends  of  the  pressure  amplitu¬ 
des,  the  disc  force  amplitude  curves  are  flat  for  on  larger  than  about  12,  then 
there  are  noticeable  increases  as  the  cavitation  number  is  decreased.  This 
characteristic  is  linked  to  the  appearance  and  growth  of  blade  sheet  cavities. 

When  comparing  the  first  blade  rate  harmonic  of  disc  force  for  the  two  dif¬ 
ferent  tip  clearance  cases,  the  force  amplitude  for  az/D  =  0.26  is  about  20  per¬ 
cent  larger  than  the  force  amplitude  for  az/D  =  0.41.  This  relative  difference 
of  force  amplitude  is  considerably  less  than  that  for  the  pressure  amplitude  at 
the  disc  center  for  the  same  tip  clearance  ratios.  For  instance,  from  data 
displayed  in  Figures  11  and  12,  at  <?n  =  8.36,  the  blade  rate  pressure  amplitude 
for  the  smaller  tip  clearance  ratio  is  170  percent  larger  than  that  for  az/D  = 
0.41.  Thus,  the  variation  of  pressure  pulse  amplitude  versus  cavitation  number 
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measured  at  one  surface  point  location  is  not  necessarily  indicative  of  the 
magnitude  of  force  amplitude  variation  acting  on  a  surrounding  reference  area 
for  this  propeller  and  wake. 

Simultaneously  with  the  disc  force  measurements,  pressure  amplitude  and 
phase  angle  measurements  were  obtained  at  11  locations  on  the  disc.  For  the 
first  three  blade  rate  harmonics  components,  the  measured  pressure  values  were 
integrated  over  the  surface  of  the  flat  disc  to  determine  a  net  normal  force 
estimate  on  the  disc  area.  Proper  account  was  made  for  the  variation  of  phase 
angle  for  each  transducer  location,  so  that  the  presence  of  transverse  asym¬ 
metries  of  pressure  amplitude  distributions  should  not  affect  the  accuracy  of 
the  integrated  pressure  results.  The  results  provide  interesting  comparisons 
with  the  measured  force  amplitudes  at  various  cavitation  numbers  at  the  two  dif¬ 
ferent  tip  clearances. 

Figures  27  and  28  compare  the  blade  rate  harmonic  force  coefficient  values 
of  measured  force  and  integrated  pressure  force  for  the  tip  clearance  ratios 
az/D  ■  0.26  and  0.41,  respectively.  The  comparisons  of  twice  blade  rate  harmonic 
force  coefficients  for  the  two  tip  clearances  are  given  in  Figures  29  and  30; 
and  the  three  times  blade  rate  amplitude  coefficients  are  displayed  in  Figures 
31  and  32,  for  az/D  ■  0.26  and  0.41,  respectively.  In  each  of  the  comparisons, 
for  all  tne  harmonic  components,  the  integrated  pressure  force  amplitudes  exceed 
the  actual  measured  disc  force  amplitudes.  The  '•'xscrepancy  is  larger  for  the 
blade  rate  harmonic  components  than  for  the  higher  harmonics.  This  serves  as 
another  indication  that  a  single  point  surface  pressure  amplitude  is  not 
necessaiily  an  accurate  representation  of  the  net  force  acting  over  a  reference 
area  surrounding  the  point. 

Comparison  with  Previous  Results 

Experimental  results  for  the  pressure  pulse  excitation  were  obtained  pre¬ 
viously  for  propeller  model  4677  in  the  large  cavitation  tunnel  at  SSPA,^  where 
a  complete  scale  model  of  the  AO-177  hull  with  appendages  was  installed  in  order 
to  produce  the  desired  model  wake.  Unsteady  pressure  amplitude  measurements  were 
recorded  at  several  locations  on  the  hull  surface  near  the  propeller.  The 
nominal  wake  in  the  water  tunnel  flow  behind  the  AO-177  model  was  not  measured, 
but  was  assumed  to  match  sufficiently  well  the  nominal  wake  measured  in  a  towing 
basin  behind  a  model  having  the  same  scale. 
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Figure  33  shows  a  comparison  of  the  longitudinal  distribution  of  the  blade 
rate  pressure  amplitude  coefficients  along  the  flat  surface'  from  the  present 
24-inch  water  tunnel  results  and  along  the  hull  surface  from  the  complete  model 
hull  tests  in  the  SSPA  tunnel.  Although  this  is  not  completely  consistent 
because  of  the  difference  between  the  flat  plate  boundary  and  the  actual  hull 
shape  and  the  mismatch  of  pressure  gauge  locations,  it  does  provide  an  approxi¬ 
mate  comparison  of  the  general  trend  of  the  longitudinal  distribution  of 
pressure  pulse  amplitudes.  Measurements  taken  on  the  flat  plate  boundary  posi¬ 
tioned  with  tip  clearance  ratios  az/D  *  0.26  and  0.41  are  shown  in  Figure  33 
together  with  measurements  obtained  with  the  model  hull  shape  with  the  vertical 
tip  clearance  of  az/D  -  0.292,  all  with  the  same  cavitation  number  an  *  8.36. 

The  values  from  the  actual  hull  shape  with  az/D  *  0.292  fall  generally  between 
the  values  obtained  with  the  flat  plate  at  the  two  different  tip  clearances. 

This  indicates  that  the  measurements  on  a  flat  plate  represent  a  good  order  of 
magnitude  approximation  for  the  measurements  along  the  hull.  However,  the  shape 
of  the  curve  of  the  (Kp)z  values  versus  xp^/D  is  not  the  same  for  the  flat  plate 
cases  and  the  hull,  and  the  details  of  differences  are  attributable  to  the  dif¬ 
ference  in  wake  and  boundary  geometry  (hull  shape). 

The  variations  wit*  respect  to  cavitation  number  of  the  blade  rate  com¬ 
ponent  of  pressure  amplitudes  at  or  near  the  tip  plane  are  shown  in  Figure  34. 
The  values  measured  on  the  complete  hull  shape  model  were  taken  over  the  pro¬ 
peller  tip  and  the  values  measured  on  the  flat  plate  boundary  were  obtained  in 
line  with  the  propeller  center.  The  curves  compare  favorably  in  an  approximate 
way,  considering  the  differences  in  boundary  shape  and  tip  clearances. 

DISCUSSION 

Useful  measurements  are  described  for  three  blade  rate  harmonics  of  the 
unsteady  excitating  force  and  pressures  on  a  reference-area  disc  mounted  in  a 
flat  boundary  above  a  cavitating  propeller  operating  in  a  screen-generated  wake 
in  the  24-inch  variable  pressure  water  tunnel  at  DTNSRDC.  The  inflow  velocity 
distribution  was  intended  to  simulate  the  axial  velocity  pattern  of  a  very 
severe  single-screw  ship  type  of  wake,  similar  to  that  of  the  AO-177.  Results 
of  measurements  are  presented  for  a  range  of  cavitation  numbers  covering 
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non-cavitating  and  cavitating  conditions  and  for  two  different  tip  clearances. 
Comparisons  are  provided  between  the  measured  disc  force  and  the  force  obtained 
by  integrating  the  pressures  over  the  disc  area.  The  integrated  pressure  force 
values  are  always  somewhat  larger  than  the  measured  force  values,  principally 
because  of  the  lack  of  sufficient  detail  in  defining  the  variation  of  phase 
angles  for  the  pressures  over  the  disc. 

Measured  results  from  the  complete  array  of  pressure  transducers  over  the 
flat  plate  boundary  show  the  longitudinal  and  transverse  centerline  distribu¬ 
tions  of  pressure  pulse  amplitudes  in  the  region  around  the  propeller.  In  general 
these  show  that  the  largest  pressure  excitation  amplitudes  occur  near  the  pro¬ 
peller  centar. 

All  the  computed  harmonics  of  unsteady  disc  force  and  pressures  show  a 
noticeable  increase  with  decreasing  cavitation  number,  particularly  below  the 
onset  of  visible  blade  sheet  cavitation  at  around  an  «  12.  Unsteady  force  and 
pressure  amplitudes  are  largest  for  the  smaller  of  the  two  tip  clearance  values. 

A  reciprocity  experiment  was  performed  to  measure  the  pressure-to-accelera- 
tion  transfer  function  (p/a),  needed  to  estimate  the  time  rate  of  change  of 
cavity  volume.  The  transfer  function  was  found  to  have  a  noticeable  variability 
at  frequencies  below  about  300  Hz,  and  this  was  judged  to  represent  an  uncer¬ 
tainty  in  the  capability  to  approximate  the  value  of  cavity  volume  velocity.  As 
seen  in  Appendix  B,  the  estimated  blade  rate  cavity  volume  velocity  inferred 
from  reciprocity  measurements  and  the  calculated  values  agree  in  order  of  magni¬ 
tude  and  trend  with  on-variation.  The  experimental  values  are  higher  than  the 
theory,  and  this  could  he  due  to  the  role  of  unsteady  tip  vortex  cavitation, 
which  is  not  accounted  for.  It  is  not  possible  at  this  time  to  sort  out  the 
magnitude  of  this  effect  on  the  discrepancy,  which  includes  the  combined 
uncertainty  of  the  experimental  calibrations  and  the  sensitivity  of  computed 
cavity  volume  to  the  propeller  inflow. 

Comparisons  between  the  present  measured  unsteady  pressure  amplitudes  pro¬ 
duced  with  the  AO-177  model  propeller  on  a  flat  plate  boundary  and  those  obtained 
previously  on  a  complete  scale  model  hull  at  SSPA  show  reasonable  correlation 
with  respect  to  the  effect  of  tip  clearance  at  comparable  cavitation  numbers. 

The  differences  between  the  results  in  terms  of  spatial  distribution  and 
variation  with  cavitation  number  are  associated  with  the  detailed  differences  in 
wake  and  boundary  geometry. 
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It  should  be  noted  that  with  the  capability  developed  here,  propeller 
excitation  experiments  could  be  carried  out  with  partial  stern  or  dummy  model 
representations  of  the  near-propeller  hull  geometry  that  include  charac¬ 
terizations  of  the  reference-area  unsteady  vertical  surface  force  as  well  as 
surface  pressures.  The  disc  force  measuring  system  described  here  could  be  used 
to  measure  the  vertical  surface  force  amplitudes  on  a  shaped  element  imbedded  in 
the  surface  of  a  foreshortened  hull  model.  Pressure  gauges  could  be  installed 
flush  with  the  shaped  surface  as  usual.  The  use  of  a  partial  body  would  provide 
some  help  with  the  simulation  of  a  desired  wake  configuration,  and  it  could  be 
made  less  subject  to  low  tunnel  speed  restrictions  Imposed  by  the  use  of  wake 
screens  a? one. 
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EXIT 

NOZZLE 


Figure  2  -  Schematic  of  Force  Measuring  Device  Installed  in 
the  24-Inch  Water  Tunnel  ODen  Jet  Test  Section 
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Figure  4  -  Calibrated  Dynamic  Response  Curve  for  Force  Measuring  Device 


PRESSURE-TO-ACCELERATION  TRANSFER  FUNCTION 


Figure  5  -  Reciprocity  Calibration  for  the  Driven  Force  Disc: 

Pressure-to-Acceleration  Transfer  Function  versus  Frequency 
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Figure  6  -  Drawings  of  Propeller  Model  4677 
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0  Towing  Basin  Wake,  Model  Hull  5326-1 

□  Screen-Wake  in  Water  Tunnel,  a  /D  -  0.26 

z 

A  Screen-Wake  in  Water  Tunnel,  a  / D  ■  0.41 
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Figure  8  -  Circumferential  Distribution  of  Axial  Wake  Velocity 
Ratio,  Comparison  of  Measurements  in  Water  Tunnel 
and  Towing  Basin  at  r/R  =  0.557 


VELOCITY  RATIO 


Figure  9  -  Circumferential  Distribution  of  Axial  Wake  Velocity 
Ratio,  Comparison  of  Measurements  in  Water  Tunnel 
and  Towing  Basin  at  r/R  =  0.774 
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Figure  11  -  Three  Blade  Rate  Harmonics  of  Measured  Pressure  Amplitude 
Coefficients  Taken  at  the  Disc  Center  with  a  /D  =  0.26 
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CAVITATION  NUMBER,  0 


Figure  12  -  Three  Blade  Rate  Harmonics  of  Measured  Pressure  Amplitude 

Coefficients  Taken  at  the  Disc  Center  with  a  /D  =  0.41 
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BLADE  RATE  PRESSURE  AMPLITUDE  COEFFICIENT 


a  /D  =  0.26 
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Figure  13  - 


Blade  Rate  Harmonic  of  Measured  Pressure  Amplitude 

Coefficient  versus  Longitudinal  Position  for  Five 

Cavitation  Numbers,  with  a  /D  =  0.26 
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THRICE  BLADE  RATE  PRESSURE  AMPLITUDE  COEFFICIENT 


Figure  13  -  Thrice  Blade  Rate  Harmonic  of  Measured  Pressure  Amplitude 

Coefficient  versus  Longitudinal  Position  for  Five  Cavitation 

numbers,  with  a  /D  =  0.26 
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TWICE  BLADE  RATE  PRESSURE  AMPLITUDE  COEFFICIENT 
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Figure  17  -  Twice  Blade  Harmonic  of 'Measured  Pressure  Amplitude 

Coefficient  veibuf.  Longitudinal  Position  for  Five 

Cavitation  Number b,  with  a  /D  =  0.41 
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TRANSVERSE,  y/D 


Figure  20  -  Twice  Blade  Rate  Harmonic  of  Meisured  Pressure  Amplitude 
Coefficient  versus  Transverse  Position  for  Several 
Cavitation  Numbers,  with  a  /D  =  0.26 
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Figure  22  -  Blade  Rate  Harmonic  of  Measured  Pressure  Amplitude 
Coefficient  versus  Transverse  Position  for  Five 
Cavitation  Numbers,  with  a  /D  =  0.41 
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Figure  26  -  Three  Blade  Rate  Harmonics  of  Measured 

Disc  Force  Coefficient  with  a  /D  =  0.41 
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BLADE  RATE  SURFACE  FORCE  COEFFICIENT 
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Figure  28  - 


CAVITATION  NUMBER,  0 


Blade  Rate  Harmonic  of  Disc  Force  Coefficient  versus 

Cavitation  Number,  Comparison  Between  Measured  Force 

and  Integrated  Pressure,  with  a  /D  =  0.41 
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Figure  33  -  Comparison  of  Measured  Blade  Rate  Pressure 
Amplitude  Longitudinal  Distribution 
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Figure  35  -  Schematic  of  the  Measurement  of  Cavit.y  Volume  Velocity  by  Reciprocal  Method 


TABLE  1  -  SUMMARY  OF  INSTRUMENTATION  CHARACTERISTICS 


Force 

Gauges 

Accelerometer 

Pressure  Gauges 

Make /Type 

Wilcoxon 

Model  L10 
Piezoelectric 

Edevco 

Model  2217E 

CEC 

Type  4-312 

Strain  gauged 
diaphragm 

Kulite 
XTM-190 
Strain  gauged 

Range 

down  to  10“^  lb 

to  250g 

to  25  psi 
(172  kPa) 

to  10  psi 
(69  kPa) 

Signal 

Condition¬ 

ing 

Gilmore 

Zero  Drive 

Gilmore 

Zero  Drive 
Model  N400 

Vi shay 

Vishay 

Filter 

1000  Hz 

5  kHz 

1000 

Hz 
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V2  is  nade  at  point  2  in  the  sane  direction  as  the  calibration  drive  force, 
with  the  cavltating  propeller  running  as  the  system  excitation.  The  volume  velo¬ 
city  of  the  cavity  at  point  1  is  thus  obtained  as 


(A. A) 


For  harmonic  oscillations  of  the  system,  the  acceleration  at  point  2  in 
the  same  direction  as  the  velocity  V2  is 


iwv. 


(A. 5) 


Thus,  another  version  of  the  basic  reciprocity  relationship  Is  expressed  as 


iw  p, 


(A. 6) 


Note  that  the  variables  measured  for  the  calibration  and  final  measurement  phases 
can  be  rearranged  for  convenience.  For  example,  suppose  that  during  the  calibra¬ 
tion  phase  the  acceleration  &2  at  the  drive  point  is  measured  as  well  the  cor¬ 
responding  fluid  pressure  pj.  In  this  case,  the  appropriate  transfer  function 
is  the  ratio  (pi/a2).  Then  the  run  measurement  of  F2  (in  the  same  direction  as 
the  calibrating  drive  acceleration)  leads  to  the  determination  of  the  cavity 
volume  velocity  from  Equation  (A. 6).  In  either  of  che  forms,  the  calibrated 
transfer  function  is  a  complex  function  of  the  real-valued  frequency  <*>.  The 
other  measured  quantities  involved  must  also  be  treated  as  complex. 

The  basic  form  of  the  simplified  reciprocity  relationship  of  Equation  (A .4) 
or  (A. 6)  can  be  derived  in  a  formal  way  using  the  second  Green  theorem  following 
concepts  discussed  by  Chertock.^ 
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APPENDIX  A 


RECIPROCITY  CONCEPT  AND  APPLICATION 

The  reciprocity  principle  states  that  for  a  linear  system,  there  is  a 
reciprocal  relationship  between  the  force  applied  and  the  system  response  at  any 
two  points.  That  is,  the  system  vibration  velocity  response  v,  say,  at  point  a 
due  to  a  force  applied  at  point  2  is  the  same  as  the  response  at  point  2  due  to  a 
force  at  point  1.  Thus, 


(A.l) 


or 


=  V  TT 
2  2 


(A. 2) 


For  the  assumed  linear  system,  the  resulting  flow  of  power  monitored  at  the  check 

points  can  be  expressed  in  terms  of  other  conjugate  variables,  for  example, 

0 

volume  velocity  -V-  times  fluid  acoustic  pressure  p. 


Then 


(A. 3) 


Figure  35  provides  a  definition  sketch  and  outline  for  the  experiments  needed  to 
make  the  reciprocity  calibration  and  to  provide  one  of  several  interpretations 
of  the  expression  given  above. 

For  the  calibration  phase,  a  mechanical  shaker  can  be  used  to  provide  an 
input  force  at  point  2  on  a  body  surface  element  (or  within  the  hull  of  a 
complete  body).  The  fluid  pressure  is  measured  simultaneously,  at  point  1  in 
the  propeller  plane  where  the  center  of  the  unsteady  cavitation  volume  is 
expected  to  occur.  A  transfer  function  is  formed  from  the  ratio  (F2 /pi )  and  this 
remains  as  a  property  of  the  system.  Then  the  measurement  of  vibration  velocity 


TABLE  2  -  GEOMETRIC  CHARACTERISTICS  OP  AO-177  CLASS  AUXILIARY 
OILER  HIGHLY  SKEWED  PROPELLER  DESIGN 


blade  chord 
diameter 

pitch 

diameter 

skew 

angle 

thickness 

chord 

camber 

chord 

rake 

diameter 

r/R 

c/D 

P/D 

9s<deg) 

/  c 

fm/c 

total 

0.2 

0.207 

1.125 

0.0 

0.20 

0.0490 

0.0 

0.3 

0.2456 

1.223 

2.2 

0.1625 

0.0444 

0.0058 

0.4 

0.2722 

1.288 

7.1 

0.1325 

0.0367 

0.021 

0.5 

0.2817 

1.318 

13.1 

0.1080 

0.0314 

0.041 

0.6 

0.2684 

1.309 

20.0 

0.0880 

0.030 

0.065 

0.7 

0.2320 

1.250 

27.7 

0.0715 

0.0295 

0.0913 

0.8 

0.1815 

1.140 

34.5 

0.0590 

0.0281 

0.109 

0.9 

0.1180 

0.970 

40.3 

0.0500 

0.0263 

0.1168 

1.0 

0.0 

0.722 

45.0 

0.045 

0.0240 

0.1148 

Number  of  blades 

Section  meanline 

Section  thickness  distribution 

Expanded  area  ratio 

Mean  width  ratio 

Blade  thickness  fraction 

Projected  skew  at  tip 

Model  diameter 

Full  scale  diameter 


7 

a*0.8 

NACA  66  Modified 
0.77 
0.216 
0.062 
45  deg 

9.812  in. (24. 92  cm) 
21  ft  (6.401  m) 


APPENDIX  B 

CAVITY  VOLUME  VELOCITY  MEASUREMENTS  AND  COMPARISONS 


Disc  force  measurements  obtained  with  excitation  from  the  operating  pro¬ 
peller,  together  with  the  calibrated  transfer  function  (p/a)  shown  in  Figure  5 
have  been  used  to  make  estimates  of  the  cavity  volume  velocity  harmonic  from  the 
expression  in  Equation  (2).  The  unsteady  disc  force  due  to  the  effect  of  cavity 
volume  variation  must  be  isolated  from  the  measured  total  force  variation 
by  vectorially  subtracting  away  the  noncavitating  or  reference  force,  taking 
proper  account  of  amplitude  and  phase.  At  the  blade  rate  frequency  of  98  Hz  for 
the  seven-b laded  propeller,  an  average  value  of  the  pressure-to-acceleration 
transfer  function  was  determined  to  oe 

(p/a)  -  0.6747  slug/ft2  (105.99  kg/m2) 

This  is  for  the  case  of  a2/D  *  0.26,  for  which  the  calibrated  transfer  function 
of  Figure  5  applies. 

The  blade  rate  disc  force  amplitude  coefficient  and  phase  angle  for  the 
noncavitating  propeller  at  az/D  -  0.26  were  found  to  be  (103)Kpnc  -  0.1622  and 
4>nc  "  Then  Table  B.l  provides  the  amplitude  and  phase  data  for  the  disc 

force  characteristics  used  to  obtain  the  inferred  blade  rate  harmonic  of  cavity 
volume  velocity  given  in  nondimens ional  form. 

Measured  and  computed  blade  rate  results  for  volume  velocity  coefficient  versus 
cavitation  number  are  compared  in  Figure  B.l,  where  calculated  values  from  two 
versions  of  the  propeller  analysis  scheme  PUF-3  are  displayed.  The  original  M.I.T. 
cavitating  propeller  analysis  program  developed  by  Lee^®  is  designated  here  as 
PUF-3(0LD).  Impr  >vements  and  non-linear  leading  edge  corrections  to  the  linear 
cavitating  blade  !  ection  analysis  have  been  incorporated  into  the  original  com¬ 
putation  scheme,  a?  outlined  In  the  paper  by  Kerwin,  Kinnas,  Wilson,  and  McHugh. ^ 9 
Results  of  the  impioved  analysis  are  denoted  by  PUF-3(NEW).  It  can  be  seen  that  the 
revised  theory  predicts  smaller  sheet  cavity  volumes  than  does  the  original  calcula¬ 
tion.  Thi3  coincid>  s  with  the  prediction  of  shorter  cavity  lengths  as  well.  The 
measured  values  of  cavity  volume  also  Include  the  effect  of  substantial  fluctuating 
tip  vortex  cavitation,  a  contribution  not  included  in  the  analytical  result. 
Considering  the  various  approximations  and  experimental  uncertainties  involved,  the 
comparison  shown  in  Figure  B.l  is  reasonable. 
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TABLE  B.l  -  ESTIMATE  OF  BLADE  RATE  CAVITY  VOLUME  VELOCITY 
USING  RECIPROCAL  MEASUREMENTS 


0 

n 

Total 

Coeff 

Oo’JKft 

Force 

Phase 

$x(deg) 

Cavltatlng  Force 
Coeff  Phase 

(103)KFc  M<*eg> 

Cavity  Volume 
Velocity  Coeff. 

<*cV"r3 

10.03 

2.052 

-16.37 

1.899 

17.91 

0.00081 

9.47 

3.538 

40.6 

3.414 

42.3 

0.00146 

8.36 

5.22 

-4.3 

5.06 

4.46 

0.00216 

i 


j 
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Measured 

Calculated,  PUF-3(0LD) 

Calculated,  PUF-3(NEW)  with  L.E. Correct ions 


Figure  B.l  -  Comparison  of  Measured  and  Calculated  Blade  Rate 
Cavity  Volume  Velocity  versus  Cavitation  Number 
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